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presence of a fully static nuclear spin system, with no ensemble averaging over initial conditions. 
In contrast to an earlier semiclassical calculation, which neglects the exchange interaction, we find 
that the singlet-triplet correlator shows a long-time saturation value that differs from 1/2, even 
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by a universal (parameter independent) value of 3n/4 at long times. The oscillation frequency 
and time-dependent phase shift of the singlet-triplet correlator can be used to perform a precision 
measurement of the exchange interaction and Overhauser field fluctuations in an experimentally 
accessible system. We also address the effect of orbital dephasing on singlet-triplet decoherence, 
and find that there is an optimal operating point where orbital dephasing becomes negligible. 
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I. INTRODUCTION 

Decoherence due to the coupling of a qubit to its en- 
vironment is widely regarded as the major obstacle to 
quantum computing and quantum information process- 
ing in solid-state systems. Electron spins confined in 
semiconductor quantum dotsi couple to their environ- 
ments primarily through the spin-orbit interaction and 
hypcrfinc interaction with nuclear spins in the surround- 
ing lattice To reach the next step in coherent electron 
spin state manipulation, the strongest decoherence ef- 
fects in this system must be understood and reduced, if 
possible. 

The effects of spin-orbit interaction are reduced in con- 
fined quantum dots at low temperatures^ Indeed, re- 
cent experiments give longitudinal relaxation times Ti for 
quantum-dot-confined electrons that reach T\ ss 20 ms 5 
in self-assembled dots and Ti k, 0.85 ms in gated dots&, 
in agreement with theoryii These times suggest that the 
spin-orbit interaction is a relatively weak source of de- 
coherence in these structures since theory predicts that 
the transverse spin decay time T 2 due to spin-orbit in- 
teraction alone (neglecting other sources of decoherence) 
would be given by T 2 = 2Ti i Other strategies for reduc- 
ing the effects of spin-orbit interaction may include using 
hole (instead of electron) spin, where a recent study has 
found that T 2 = 2Ti also applies, and the hole spin re- 
laxation time can be made even longer than that for the 
electron spin£ 

Unlike the spin-orbit interaction, the hyperfine inter- 
action of a single electron spin with a random nuclear 
spin environment can lead to pure dephasing, giving a 
transverse spin decay time on the order of 5 nsj^ilfiiAi 
six orders of magnitude shorter than the measured lon- 
gitudinal decay times T±. To minimize errors during 
qubit gating operations in these proposed devices, this 
decay must be fully understood. The hyperfine interac- 



tion in a single quantum dot is described by a Hamil- 
tonian H = h ■ S, where S is the electron spin opera- 
tor and h is a collective quantum nuclear spin operator, 
which we will refer to as the "Overhauser operator" . A 
common assumption in the literature is to replace the 
Overhauser operator by a classical effective magnetic field 

h -> B^VVVWVVVVWV 2 SWCB & daSSiCal 

magnetic field only induces precession (not decoherence) , 
the classical-field picture necessitates an ensemble of nu- 
clear spin configurations to induce decay of the electron 
spin expectation valueA^S For experiments performed on 
a large bulk sample of electron spins, or experiments per- 
formed over timescales that are longer than the typical 
timescale for variation of Bat, the source of the ensemble 
averaging is clear. However, one conclusion of this model 
is that single-electron-spin experiments performed over a 
timescale shorter than the nuclear spin correlation time 
should show no decay. This conclusion is contradicted 
by numericalSi^ and analyticalii^ results, which show 
that the quantum nature of the Overhauser operator can 
lead to rapid decay of a single electron spin, even for 
a fully static nuclear spin system. This rapid decay is, 
however, reversible with a standard Hahn spin-echo se- 
quence in an applied magnetic field and the timescale of 
the decay can be increased by squeezing the nuclear spin 
stated 

Another potential solution to the hyperfine decoher- 
ence problem is to polarize the nuclear spins. Polarizing 
the nuclear spin system in zero applied magnetic field 
reduces the longitudinal spin-flip probability by the fac- 
tor l/p 2 N , where p is the nuclear spin polarization and 
N is the number of nuclear spins within the quantum 
dot£ii The effect on the transverse components of elec- 
tron spin is different. Unless the nuclear spin state is 
squeezed or a spin-echo sequence is performed, the trans- 
verse components of electron spin will decay to zero in 
a time t c ss 5 ns in a typical GaAs quantum dot. Po- 
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larizing the nuclear spin system increases t c by reducing 
the phase-space available for fluctuations in the Over- 
hauser operator, resulting in t c w 5 ns/ y/1 — p 2 m^ Re- 
cent experiments show that the nuclear spin system can 
be polarized by as much as 60%>i& However, to achieve an 
order-of-magnitude increase in t c , the polarization degree 
would have to be on the order of 99%^ for which more 
ambitious polarization schemes have been proposed^ 

If electron spins in quantum dots are to be used as 
quantum information processors, the two-electron states 
of double quantum dots must also be coherent dur- 
ing rapid two-qubit switching times4& Measurements 
of singlet-triplet relaxation times t$T in vertical dou- 
ble dots (tsr ~ 200 fis) gated lateral double dots 
(t S T ~ 70/is)r& and single dots (t ST ~ 2.58ms)Si sug- 
gest that these states may be very long-lived. Recent ex- 
periments have now probed the dccohcrcnce time of such 
states, which is believed to be limited by the hyperfine 
interaction with surrounding nuclear spins^S The dra- 
matic effect of the hyperfine interaction on two-electron 
states in a double quantum dot has previously been il- 
lustrated in experiments that show slow time-dependent 
current oscillations in transport current through a double 
dot in the spin blockade regime^ 

It may be possible to circumvent some of the complica- 
tions associated with single-spin decoherence by consid- 
ering an encoded qubit, composed of the two-dimensional 
subspace of states with total z-projection of spin equal 
to zero for two electrons in a double quantum dotiiS One 
potential advantage of such a setup is that it may be 
possible to reduce the strength of hyperfine coupling to 
the encoded state space for a symmetric double-dot (see 
Appendix^). A potential disadvantage of this scheme is 
that coupling to the orbital (charge) degree of freedom 
can then lead to additional decoherence, but we find that 
orbital dephasing can be made negligible under appropri- 
ate conditions (see Sec. II VL To achieve control of the 
singlet-triplet subspace, however, the decoherence pro- 
cess for the two-electron system should be understood in 
detail. 

In this paper we give a fully quantum mechanical so- 
lution for the spin dynamics of a two-electron system 
coupled to a nuclear-spin environment via the hyperfine 
interaction in a double quantum dot. Although we focus 
our attention here on quantum dots, decoherence due to 
a spin bath is also an important problem for, e.g., pro- 
posals to use molecular magnets for quantum information 
processing j 3 i 31 i 32 i 33 In fact, the problem of a pair of elec- 
trons interacting with a bath of nuclear spins via the con- 
tact hyperfine interaction has been addressed long ago to 
describe spin-dependent reaction rates in radicals>i2iM A 
semiclassical theory has been developed^ in which elec- 
tron spins in radicals experience a randomly oriented ef- 
fective classical magnetic field due to the contact hyper- 
fine interaction between electron and nuclear spins. In 
this semiclassical theory, random hopping events of the 
electrons were envisioned to induce a randomly fluctuat- 
ing local magnetic field at the site of the electron spin, 



resulting in decay of a singlet-triplet correlator. Here, 
we solve a different problem. Ensemble averaging over 
nuclear spin configurations is natural for a large sample 
of ~ 10 23 radicals. In contrast, we consider the coher- 
ent dynamics of two-electron spin states within a single 
double quantum dot. More importantly, the Heisenbcrg 
exch ang e interaction, which was found to be negligible in 
Ref. Il2l can be any value (large or small) in our system 
of interest. We find that a nonzero exchange interaction 
can lead to a drastic change in the form and timescale 
of decoherence. Moreover, this paper is of direct rele- 
vance to very recent experiments22iSi2ai related to such 
double-dot systems. 

The rest of this paper is organized as follows. In Sec. 
[n]we solve the problem for electron spin dynamics in the 
subspace of total spin z-component S z — with an ex- 
act solution for the projected effective Hamiltonian. In 
Sec. IIIII we show that a perturbative solution is possible 
for electron spin dynamics in the subspace of singlet and 
S z = +1 triplet states. Sec. El contains a discussion of 
the contributions to singlet-triplet decoherence from or- 
bital dephasing. In Sec. 0we review our most important 
results. Technical details are given in Appendixes [S] to 

n 

II. DYNAMICS IN THE S z = SUBSPACE 

We consider two electrons confined to a double quan- 
tum dot, of the type considered, for example, in Rcfs. 
I20ll2lll3l Each electron spin experiences a Zeeman split- 
ting e 2 = gfiBB due to an applied magnetic field B = 
(0,0, B), B > 0, defining the spin quantization axis z, 
which can be along or perpendicular to the quantum dot 
axis. In addition, each electron interacts with an inde- 
pendent quantum nuclear field h;, I = 1,2, due to the 
contact hyperfine interaction with surrounding nuclear 
spins. The nuclear field experienced by an electron in 
orbital state I is h; = ^k^-k, where is the nuclear 
spin operator for a nucleus of total spin / at lattice site fc, 
and the hyperfine coupling constants are given by A\ = 

vA \ ipo( T k)\ 2 , with v the volume of a unit cell containing 
one nuclear spin, A characterizes the hyperfine coupling 
strength, and t/j l (rk) is the single-particle envelope wave- 
function for orbital state I, evaluated at site k. This prob- 
lem simplifies considerably in a moderately large mag- 
netic field (B > max{((5h) rms /g/i B , (h) rms /g£t B }, where 

(O) = O 2 1%/jj) 1 ^ 2 is the root-mean-square expec- 
tation value of the operator O with respect to the nuclear 
spin state Sh = ^ (hi — h^), and h = | (hi + h.2)). 
In a typical unpolarizcd GaAs quantum dot, this con- 
dition is B 1 Aj \J~N g^LB ~ 10 mT (see Appendix 
|A"|) . For this estimate, we have used IA/g/iB « 5T, 
based on a sum over all three nuclear spin isotopes (all 
three hyperfine coupling constants) present in GaAs^i 
and N « 10 5 nuclei within each quantum dot. In this 
section, we also require B J/g/is, where J is the 
Heisenberg exchange coupling between the two electron 
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spins. For dcfiniteness we take J > 0, but all results are 
valid for either sign of J, with J replaced by its abso- 
lute value. In the above limits, the electron Zeeman en- 
ergy dominates all other energy scales and the relevant 
spin Hamiltonian becomes block-diagonal, with blocks 
labeled by the total spin projection along the magnetic 
field S z (see Appendix [BJ). In the subspace of S z = we 
write the projected two-electron spin Hamiltonian in the 
subspace of singlet and S z — triplet states (\S) , |T )) 
to zeroth order in the inverse Zeeman splitting l/e 2 as 
H = jS-S + 5h z 5S z , where S = S! + S 2 is the total spin 
operator in the double dot and SS = Si — S2 is the spin 
difference operator. In terms of the vector of Pauli matri- 
cesr = (t*,t^,t z ): \S) -> \t z = -1) , \T ) -> \t z = +1) 
Hq can be rewritten as: 



Hq ='i(l + T z ) + 6h z T x 



(1) 



Diagonalizing this two-dimensional Hamiltonian gives 
eigenvalues and eigenvectors 



\Et) 



J -± l - 

2 2 



J 2 + 4 (6h*) 
5h z n \S)+E±\T ( 



(g) \n) 



(SKY 



(2) 
(3) 



where \n) is an eigenstate of the operator 5h z with eigen- 
value Sh z n . Since the eigenstates \E^) are simultaneous 
cigenstatcs of the operator 5h z , we note that there will 
be no dynamics induced in the nuclear system under the 
Hamiltonian Hq. In other words, the nuclear system re- 
mains static under the influence of Hq alone, and there 
is consequently no back action on the electron spin due 
to nuclear dynamics. 

We fix the electron system in the singlet state \S) at 
time t = 0: 



w=o)) = \s)®\i>i); m=J2' 



(4) 



where a n is an arbitrary set of (normalized) coefficients 
(En l a «| 2 = !)• The initial nuclear spin state \ipj) is, in 
general, not an eigenstate \n). The probability to find 
the electron spins in the state \Tq) at t > is then given 
by the correlation function (setting H = 1): 



C To (t) = yPi(n) 



®{T \ 



-,-iHnt 



(5) 



where pi(n) 



gives the diagonal matrix elements of 



the nuclear-spin density operator, which describes a pure 
(not mixed) state of the nuclear system: pj = \ipj) = 

E„ PliP) \ n ) ( n \+T, n7 tn' < a rA n ') H C T„(t) is the Sum 

of a time-independent piece C„ and an interference term 



c To (t) 



°T 



(*)> 



(6) 



L n 



2 {sky 



j* + A{5h z n y 



-C n cos([E+ -En]t) 



(7) 

(8) 



Here, the overbar is defined by f(n) = EnAr( n )/( n )- 
Note that C n depends only on the exchange and Over- 
hauser field inhomogeneity Sh z n through the ratio bh z n j J . 

For a large number of nuclear spins iV > 1 in a su- 
perposition of (J/i^-eigenstates |n), we assume that pi(n) 
describes a continuous Gaussian distribution of Sh z n val- 
ues, with mean 8h^ = (for the case Sh z t ^ 0, see 

Sec. Ill All and variance a 2 = (Sh^ — 5h^ 2 = (<5/i„) 2 (i.e. 
(To = (Sh z ) rms ). The approach to a Gaussian distribution 
in the limit of large N for a sufficiently randomized nu- 
clear system is guaranteed by the central limit theorem^! 
The assumption of a continuous distribution of Sh z n pre- 
cludes any possibility of recurrence in the correlator we 
calculated A lower-bound for the Poincare recurrence 
time in this system is given by the inverse mean level 
spacing for the fully-polarized problem^: t p > N 2 /A. In 
a GaAs double quantum dot containing N ~ 10 5 nuclear 
spins, this estimate gives i rcc > 0.1s. Moreover, by per- 
forming the continuum limit, we restrict ourselves to the 
free- induction signal (without spin-echo). In fact, we re- 
mark that all decay in the correlator given by (JHJ can 
be recovered with a suitable 7r-pulsc, defined by the uni- 
tary operation U v |-E„) = \E%). This statement follows 
directly from the sequence 



-ijt 



E^) = U^- lHot U^- lH ^ \E±) 



(9) 



Thus, under the above sequence of echoes and free in- 
duction, all eigenstates are recovered up to a common 
phase factor. Only higher-order corrections to the effec- 
tive Hamiltonian Hq may induce completely irreversible 
decay. This irreversible decay could be due, for example, 
to the variation in hyperfine coupling constants, leading 
to decay on a timescale t ~ N/A, as in the case of a single 
electron spin in Refs. I9HT1I Another source of decay is 
orbital dephasing (see Sec. II VI) . 

We perform the continuum limit for the average of an 
arbitrary function f(n) according to the prescription 

$>i(n)/(n) -> f dxP^(x)f(n(x)), (10) 

n ^ 

P^(x) = ^cxpf-^-p!^), (11) 



with x = 0, a 
Using 



2 _ ^2 



2ira 



and here we take 



Shi 



a = oq. 



2x 2 



J 2 + Ax 2 ' 



C n = C{5h z n ) = C(x) = 

we evaluate Cjf(i) = Re \cf*{t) 
interference term is given by the integral 



(12) 



where the complex 



dxC(x)P tTo -o(x)e 



(13) 
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In general, the interference term given by Eq. 113|) 
will decay to zero after the singlet-triplet decoherence 
time. We note that the interference term decays even 
for a purely static nuclear spin configuration with no en- 
semble averaging performed over initial conditions, as is 
the case for an isolated electron spin»ii*SiS The total z- 
componcnt of the nuclear spins will be essentially static 
in any experiment performed over a timescale less than 
the nuclear spin diffusion time (the diffusion time is sev- 
eral seconds for nuclei surrounding donors in GaAs^I). 
We stress that the relevant timescale in the present case 
is the spin diffusion time, and not the dipolar correlation 
time, since nonsecular corrections to the dipole-dipole in- 
teraction are strongly suppressed by the nuclear Zeeman 
energy in an applied magnetic field of a few Gauss2& (as 
assumed here). Without preparation of the initial nu- 
clear state or implementation of a spin-echo technique, 
this decoherence process therefore cannot be eliminated 
with fast measurement, and in general cannot be modeled 
by a classical nuclear field moving due to slow internal 
dynamics; a classical nuclear field that does not move 
cannot induce decay. 

At times longer than the singlet-triplet decoherence 
time the interference term vanishes, leaving Ct (oo) = 
C„, which depends only on the ratio Sh^/J, and could 
therefore be used to trace-out the slow adiabatic dynam- 
ics 5h^(t) of the nuclear spins, or to measure the exchange 
coupling J when the size of the hyperfine field fluctua- 
tions is known. We evaluate Ct (oo) from 

/oo 
dxC{x)P ao , {x). (14) 
-oo 

In two limiting cases, we find the saturation value is given 
by (see Appendix Q 



C To (oo) 



~ ./ 

2 4*0 



oo > J, 

(T < J. 



(15) 



We recover the semiclassical high-magnetic-field limiU^ 
(Cjb(oo) = 1/2) only when the exchange J is much 
smaller than o~o- Furthermore, due to the average over 
b~h z n eigenstates, the approach to the semiclassical value 
of h is a slowly-varying (linear) function of the ra- 
tio J /(Jo, in spite of the fact that C n oc (J/Sh!^) 2 as 
J — > 0. In Figure we plot the correlator saturation 
value Ct (oo) as a function of the ratio (Sh z ) rms /J for 
a nuclear spin system described by a fixed eigenstate of 
Sh z (i.e. pi = \n) (n\), and for a nuclear spin system 
that describes a Gaussian distribution of Sh z eigenstates 

with variance <7q = {5h z L ) 2 = (Sh z ) 2 ms . We also show the 
asymptotic expression for (T(i> J, as given in Eq. i|15f) . 

Now we turn to the interference term C^?* if) given by 
Eq. (|13fl . which can be evaluated explicitly in several 
interesting limits. First, in the limiting case of vanishing 
exchange (J = 0), we have C(x) = \ from i|12|l . Direct 
integration of Eq. then gives 



Cjt(t)=--ex P [-- s ),to = —,J^0. (16) 
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FIG. 1: Saturation value of the triplet occupation probability 
C7t (oo) = C n vs. {Sh z ) rlas I J when the nuclear spin system 
has been squeezed into an eigenstate of 5h z (dashed line) and 
when the state of the nuclear spin system describes a Gaus- 
sian distribution of eigenvalues 8h n , with mean 5h„ = and 

variance a 2 = (Sh^) 2 (solid line). We also show the analytical 
asymptotics for co S> J, given by Eq. 11511 (dotted line) and 
the semiclassical value (Ct (oo) = 1/2) (dash-dotted line). 



For zero exchange interaction, the correlator decays 
purely as a Gaussian, with decoherence time to = ^r- ~ 

for a typical asymmetric double quantum dot (see 
Appendix 0. However, for arbitrary nonzero exchange 
interaction J / 0, wc find the asymptotic form of the 
correlator at long times is given by (see Appendix Q : 



C% (t) 



(Jt 



37T ' 

4 i 



4a VJt 3 / 2 
1 



t ^> max 



1 J 
J' 2ctq"' 4af 



(17) 
(18) 



Thus, for arbitrarily small exchange interaction J, the 
asymptotic decay law of the correlator is modified from 
the Gaussian behavior of Eq. (|16fl to a (much slower) 
power law (~ l/i 3 / 2 ). We also note that the long-time 
correlator has a universal phase shift of which is in- 
dependent of any microscopic parameters. Our calcu- 
lation therefore provides an example of interesting non- 
Markovian decay in an experimentally accessible system. 
Furthermore, the slow-down of the asymptotic decay sug- 
gests that the exchange interaction can be used to modify 
the form of decay, in addition to the decoherence time, 
through a narrowing of the distribution of eigenstates 
(see the discussion following Eq. I|2(J|) below). We have 
evaluated the full correlator CV (t) by numerical integra- 
tion of Eq. I|13fl and plotted the results in Figure |3 along 
with the analytical asymptotic forms from i|17|) . 

We now investigate the relevant singlet-triplet correla- 
tor Ct if) in the limit of large exchange J. In this case, 
we have x < cto <Si J for the typical x contributing to the 
integral in Eq. H13f) . Thus, we can expand the prcfactor 
C{x) and frequency term in the integrand: 



C{x) 



'J 2 ' 



(19) 
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FIG. 2: Decay of the correlator CT (t) evaluated by numerical 
integration of Eq. 1131 for three ratios of ao/J: ao/J = 2 
(dashed-dotted line), ao/J — 1 (dashed line), and (Jo/ J = 1/2 
(solid line). The analytical long-time asymptotic expressions 
from Eq. 11711 are shown as dotted lines. 



FIG. 3: The correlator CV (*)/Ct (°°) shows a rapid Gaus- 
sian decay when J = (solid line, from Eq. 1161 ), but has 
a much slower power-law decay ~ 1/t 3 / 2 for large exchange 
J = lOoo > 0o (dotted line, from Eq. 



\/ J 2 + 4z 2 w J + 2 



J 



(20) 



From Eq. I|20[) it is evident that the range of frequen- 
cies that contribute to the correlator is suppressed by 
(To/ J (increasing the exchange narrows the distribution of 
eigenenergies that can contribute to decay). This narrow- 
ing of the linewidth will increase the decoherence time. 
Moreover, the leading-order x 2 -dcpcndence in (|2(Jfl col- 
laborates with the Gaussian distribution of Sh z eigen- 
states to induce a power-law decay. With the approxi- 
mations in Eqs. (|19fl and H20|) . we find an expression for 
the correlator that is valid for all times in the limit of 
large exchange J by direct evaluation of the integral in 
Eq. (JT3J): 



2 cos ( Jt + | arctan 







2 x 3/4 



(21) 



(22) 



There is a new timescale (t' = J/Ao 2 /) that appears for 
large J due to dynamical narrowing; increasing the ex- 
change J results in rapid precession of the pscudospin 
r about the z-axis, which makes transverse fluctuations 
along t x due to 8h z progressively unimportant. Explic- 



we have t' 



JN/4A 2 > VN/A for J > 



CO 



itly, 
a/Vn. 

Eq. I|21|) provides a potentially useful means of extract- 
ing the relevant microscopic parameters from an experi- 
ment. J and do can be determined independent of each 
other exclusively from a measurement of the oscillation 
frequency and phase shift of C^^t). In particular, any 
loss of oscillation amplitude (visibility) due to systematic 
error in the experiment can be ignored for the purposes 
of finding erg ancl J- The loss in visibility can then be 



quantified by comparison with the amplitude expected 
from Eq. (|21|) . We illustrate the two types of decay that 
occur for large and small J in Figure [3] 



A. Inhomogeneous polarization, Sh^ ^ 

It is possible that a nonequilibrium inhomogeneous av- 
erage polarization could be generated in the nuclear spin 
system, in which case Sh^ ^ 0. Pumping of nuclear spin 
polarization occurs naturally, for example, at donor im- 
purities in GaAs during electron spin resonance (ESR), 
resulting in a shift of the ESR resonance condition, 3 ^ 
It is therefore important to investigate the effects of a 
nonzero average Overhauser field inhomogeneity on the 
decay law and timescale of the singlet-triplet correlator. 
In this subsection we generalize our previous results for 
the case bh z n ^ 0. 

We set the mean Overhauser field inhomogeneity to 
Sh^ = xq, in which case the complex singlet-triplet inter- 
ference term is given by 



C%(t) = - 



dxC(x)P ao ., Xo (x)e iWJ * +ix \ (23) 



When the mean value of the Overhauser field inho- 
mogeneity xq is much larger than the fluctuations Co 
(xq (To) i wc approximate C(x) ~ C(xq) and expand 
the frequency term V J 2 + 4x 2 = lu + ^-(x — x ) + 

^pr(x — xq) 2 +■■■■, where uio = V ' J 2 + 4xq. We retain 
only linear order in x — xq for the frequency term, which 
is strictly valid for times t < ( J 2 + 4x 2 l ) 3 / 2 /2 J 2 a 2 . This 
time estimate is found by replacing [x — xo) 2 ~ <Jq in the 
quadratic term and demanding that the quadratic term 
multiplied by time be much less than one. In this limit, 
the correlator and range of validity arc then 



cos (u>ot) . 



(24) 
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^0 

Ax a 



X » CT , t <C 



, w = ^J 2 +4a;g, (25) 
" 2 + 4x2) 3/2 



(26) 



This expression is valid for any value of the exchange J, 
up to the timescale indicated. 

In contrast with the previous result for xq = 0, from 
Eq. I|24(l we find that the long-time saturation value 
of the correlator deviates from the semiclassical result 
(C To (oo) = -C%g(0) = 1/2) by an amount that is 



quadratic in the exchange J for J <C xq ■ 



a To (oo) = c^(o) 



8 lio 
v 2 



2{J)\ J»x . 



X » cr . 



(27) 



In the limit of large exchange, J S> max (ao, Xo), we 
can once again apply the approximations given in Eqs. 
(|19|l and 12(J|) . Using these approximations in Eq. i|23|) 
and integrating then gives 



<5£(i) = - 2 (^)V(*)(l+(^) £ 2 W jexp{zJi-||(l-^))}, (28) 

^^O"^) 7 ''" = ^' J>>maX( ^ g ° M<< 2max^^) - (29) 

I 



We have found the limit on the time range of validity in 
Eq. (|29|1 using the same estimate that was used for Eqs. 
(|24I26[) . At short times, t < t' = J/4ct 2 , we expand 

£ 2 (i) « 1 + — ( jr ) and find that this function decays 

t V / 

initially as a Gaussian with timescale t ' J/Axoao- 



J 2 
J 



4x cr ' 
J 
4^ 



cos(w i), (30) 
2r 2 

= ./+^, (31) 



i <C t' Q = ^—2, J > max(xo, ctq)- (32) 



This agrees with the result in Eq. Ij24(l when J ^ xq ^> 
a . 

For sufficiently large exchange J, the expression given 
by Eq. 128|) is valid for times longer than the previous 
expression, given by Eq. (|24|l . We perform an asymptotic 
expansion of Eq. I|28(l for long times using £(t ^ t' ) ~ 
e m/4 \/%/i. This gives 



B. Reducing decoherence 

The results of this section suggest a general strategy for 
increasing the amplitude of coherent oscillations between 
the singlet \S) and triplet |T ) states, and for weakening 
the form of decay. To avoid a rapid Gaussian decay with 
a timescale t' = J/Axoao, the mean Overhauser field 
inhomogeneity should be made smaller than the fluctu- 
ations (Sh^ = xo < Co) an< 4 the exchange J should be 
made larger than xq and ao (J 3> max(xo, ao))- Explic- 
itly, the ideal condition for slow and weak (power-law) 
decay can be written as 



J 3> a > x . 



(35) 



The condition in Eq. (|35|) can be achieved equally well 
by increasing the exchange coupling J for fixed hypcrfine 
fluctuations ao or by reducing the fluctuations ao through 
state squeezing or by making the double-dot confining 
potential more symmetric (see Appendix |A"|) . 



-x 2 /2v% COS ( Jt + ) 



(33) 



AaoVJt 3 / 2 
t^>t' Q = —^, J > max(a;o, fio)- (34) 

4(T n 



As in the case of xq = 0, the long-time asymptotics of Eq. 
(|2"8|) once again give a power law ~ 1/i 3 / 2 , although the 
amplitude of the long-time decay is exponentially sup- 
pressed in the ratio x\ja^. When xo — 0, Eq. 
recovers the previous result, given in Eq. i|17|) . 



III. DYNAMICS IN THE SUBSPACE OF \S) 

AND |T+) 

We now consider the case when the Zccman energy of 
the S z = 1 triplet state approximately compensates the 
exchange (|A| <C J, where A = e z + J). In addition, 
we assume the exchange is much larger than the nuclear 
field energy scales J » max {(<5h) rms , (h) rms }. Under 
these conditions, wc consider the dynamics in a subspace 
formed by the singlet \S) — * \t z = —1) and the S z = 1 
triplet state |T + ) — > \t z = +1), governed by the Hamil- 
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tonian (to zeroth order in 1/J, sec Appendix lB|l : 
H+ = i (A + ft*) (1 +r z ) - (^"r+ + H.c.) . (36) 

Here, (5/i ± = 8h x ±i5W and r± = f (t 1 ± iY»). The \T+) 
probability at time t > is 

c T+ (t) = E^( n ) K n 'i ® < T +i e ~ 4iM i 5 > ® M 2 ■ ( 3? ) 

n,n' 

This case is essentially different from the previous one, 
since the eigenstates of H + are no longer simply prod- 
uct states of electron and nuclear spin, implying a back- 
action of the electron on the nuclear system. Never- 
theless, when (h z + A) rms » (<^h ± ) rms , we can eval- 
uate the correlator in standard time-dependent per- 
turbation theory to leading order in the term V — 
— -j= (r + 5h~ + r~5h + ). Neglecting corrections of order 

h^/A <c 1, this gives 



C T>)«^(l-cos([[n i + A]t)), (38) 

where a n = Yin' \ ( n '\ l n )| 2 ' an< ^ \ n ) m now an e ig en_ 
state of the operator h z with eigenvalue [h z ] n . To es- 
timate the size of <x n , we assume identical completely 
decoupled dots and nuclear polarization p < 1, which 
gives a n « + 1) Y^k ^k> wriere Ak is the hypcrfinc 
coupling constant to the nuclear spin at lattice site k 
(with total nuclear spin I) and the sum runs over all 
lattice sites in one of the dots. We estimate the typical 
size of a n with the replacements Ak — > -4, J2k ~ * 

which gives a n m = J 1 ^jy^ A, where N charac- 

terizes the number of nuclear spins within the dot enve- 
lope wavefunction. If we assume the nuclear spin state 
is described by a continuous Gaussian distribution of h z 
eigenstates with mean = and variance cr3_, we find 

(39) 

Thus, if we ignore any possibility for recurrence, the dis- 
tribution of h z eigenstates will lead to Gaussian decay 
of the two-electron spin state, as is the case for a single 
electroniii*Si However, as in the case of a single electron, 
this decay can be reduced or eliminated altogether by 
narrowing the distribution of h z eigenstates |n) through 
measurement (squeezing the nuclear spin state) ill We 
show these two cases (with and without squeezing of the 
nuclear state) in Figure 0] 

IV. SINGLET-TRIPLET DECOHERENCE DUE 
TO ORBITAL DEPHASING 

To this point we have neglected dephasing of the sin- 
glet \S) and triplet \Tj) (j = 0, +) states due to cou- 
pling in the orbital sector. The effective Hamiltonian 



o 











A °+ - 






■■ ■' i V .' i 





1 2 , , „ 3 4 5 

t (1/a) 

FIG. 4: Decay of the correlator Ct + (t) in two cases. A typ- 
ical unprepared initial state, where the nuclear spin system 
is in a superposition of h z eigenstates, results in a Gaussian 
decay (solid line, from Eq. (13910 . If the nuclear spin state is 
squeezed into an h z eigenstate, there is no decay, only coher- 
ent oscillations (dotted line, from Eq. 13810 . For this plot we 
have used A = 5a. 

description ignores the different character of the orbital 
states for singlet and triplet, and so it is tempting to 
assume that orbital dephasing is unimportant where the 
effective Hamiltonian is valid. However, the singlet and 
triplet do have different orbital states which can, in gen- 
eral, couple differently to the environment through the 
charge degree of freedom, and therefore acquire different 
phases. Examples of such environmental influences are 
charge fluctuators or measurement devices, such as quan- 
tum point contacts used for charge readout Here we 
briefly step away from the effective Hamiltonians derived 
in Appendix IB1 to give a physical picture of the effects of 
orbital dephasing in terms of the true double-dot wave- 
functions. We then return to the effective Hamiltonian 
picture in order to give a more general estimate of the 
effects of orbital dephasing on singlet-triplet decoherence 
for a two-electron double dot. 

We consider a double quantum dot containing a fixed 
(quantized) number of electrons N. Within the far-field 
approximation, the double-dot charge distribution cou- 
ples to the environment first through a monopolc, and 
then a dipole term. Since the charge on the double dot is 
quantized, the monopolc term gives an equal contribution 
for both the singlet and triplet wavefunctions. The lead- 
ing interaction that can distinguish singlet from triplet is 
the electric dipole term: 

Krb(i) « -Piv- E(t). (40) 

Here, pn is the electric dipole moment operator for the 
charge distribution in a double dot containing TV elec- 
trons and E(f) is a fluctuating electric field due to the 
surrounding environment, which we model by a Gaus- 
sian random process. For a double quantum dot with 
well-localized single-particle eigenstates we denote the 
charge states by |(n, m)), indicating that the double- 
dot has n electrons in dot 1 and m electrons in dot 2, 
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where n + m = N . If the double dot contains only 
a single electron (N = 1), the environment can distin- 
guish the two localized states through the difference in 
the dipole moment operator, which has the size |Api| = 
K(l,0)|pi |(1,0)) - <(0, 1)| Pl |(0, 1))| « 2\e\a, where e 
is the electron charge and 2a is the inter-dot spacing. 
When N = 2, for highly-localized states, only the states 
with double-occupancy (|(0, 2)) and |(2, 0))) contribute to 
the dipole moment. If the typical hypcrfine energy scale 
is much smaller than the detuning from resonance S of 
the |(1,1)) and |(0,2)) states (max ({Sh) rms , (h) rms ) « 
S), only the |(1,1)) singlet state (not the triplets) will 
mix with the doubly-occupied states, so the singlet and 
triplet states will be energetically distinguishable through 



|A P2 | = |<S|pa|S)| « 2\e\a P { 



(0,2) 



(2,0) 



£ 2\e\aD, 



where P( 2 ) (P(2 o) ) is the probability to find the singlet 
IS") in the |(0,2))' (|(2,0))) state and D = P, 



(0,2) 



Pi 



(2,0) 



is the double occupancy. In this discussion, we assume 
that the exchange is much larger than the hyperfine en- 
ergy scales, J > max((h) rms , (<5h) rms ), so that the sin- 
glet and triplet states are good approximates for the true 
two-electron eigenstates. 

For weak coupling to the environment, and assuming 
the environment correlation time is much less than the 
orbital dephasing time , we can apply standard tech- 
niques to determine the dephasing time for a two-level 
system described by the Bloch equations^i We find that 
the fluctuations in E(i) lead to exponential dephasing 

with the rate l/tf ] = \ |A PAr | 2 ^ oo dtE{t)E{Q), where 
the scalar E{t) is the component of E(£) along Ap^v and 
we assume lim^oo j J* dt'E(t') = 0. Assuming equiva- 
lent environments for the single-particle and two-particle 
cases, the ratio of the single-particle to two-particle de- 
phasing times is then 



,(2) 



Ap 2 



Api 



<D 2 



(41) 



The single-electron orbital dephasing rate has been mea- 
sured to be t^' ~ 1 ns^ and t^' « 400 ps^ in different 
gated double quantum dots. If the hyperfine interaction 
(which becomes important on the timescale t > 5 ns) is 
to provide the major source of decoherence in these two- 
electron structures, we therefore require 3> t)t'. This 
condition can be achieved by ensuring a small double oc- 
cupancy D <C 1 of the singlet state. When the inter-dot 
tunnel coupling t\% is much less than the detuning from 
resonance 8 (t\ 2 <C S <C U + U' , with on-site and nearest- 
neighbor charging energies U and U', respectively - see 
Appendix [5J) we find the double-occupancy of \S) in per- 
turbation theory is 



D 



< 1. 



(42) 



Even in this regime, orbital dephasing may become the 
limiting timescale for singlet-triplet decoherence after the 



removal of hyperfinc-induced decoherence by spin echo. 
A detailed analysis of the double-occupancy and its rela- 
tion to the concurrence (an entanglement meas ure) for a 
symmetric double dot can be found in Refs. lillia . 

With this physical picture in mind, we can generalize 
the above results to the case when the electrons expe- 
rience fluctuations due to any time-dependent classical 
fields. In particular, if the separation in single-particle 
energy eigenstates for N = 1 is e + Se(t), where Se(t) 
fluctuates randomly with amplitude Se, and similarly, if 
for N = 2 the singlet and triplet levels are separated by 
an exchange J + SJ(t), where 6J(t) has amplitude 5 J, 
we find 



SJ_ 
5e 



(43) 



From this expression we conclude that the optimal op- 
crating point of the double dot is where the slope of 
J vs. e vanishes, i.e., SJ/6e = 0. At this optimal 

(2) 

point, — > oo, within the approximations we have 
made. Eq. I|43|) is valid for weak coupling to the en- 
vironment (i.e. 5 J -C J and Se <C e), and when the 
environment correlation time is small compared to the 
dephasing times. If, for example, we take J w 2t\ 2 /S for 
U + U' 5 ti2 from Eq. IjBllfl and if Se corresponds 
to fluctuations in the single-particle charging energy dif- 



ference (e ~ (V g i - V g2 ) ~ 6 from Eq. {B9J), we find 

$V*i 2) ~ 4 *i2/<* 4 > in agreement with Eqs. gTJ and gg) . 
In particular, the hypcrfinc-dominated singlet-triplet de- 
coherence becomes visible when t^ 2 ' t' , i ' to, t + . 
This regime is achievable by choosing S ^S> t± 2 , but still 
J f» 2t\ 2 /5 (To, since t? is a much stronger function 
of 8 than t' , i ' . That is, the two-particle dephasing time 
scales like r?' ~ <5 4 , but the typical hyperfine- induced 
decay times scale like t' ,tQ ~ J ~ l/S. On the other 
hand, when t\ 2 ~ S, we have \SJ/6e\ ~ O(l), which gives 
*i ~ *i ' an<4 thus a very short singlet-triplet decoher- 
ence time (~ Ins), which is dominated by orbital de- 
phasing. 



V. CONCLUSIONS 

We have shown that a fully quantum mechanical solu- 
tion is possible for the dynamics of a two-electron system 
interacting with an environment of nuclear spins under 
an applied magnetic field. Our solution shows that the 
singlet-triplet correlators CT (t) and CV(i) will decay 
due to the quantum distribution of the nuclear spin sys- 
tem, even for a nuclear system that is static. We have 
found that the asymptotic behavior of Ct (i) undergoes a 
transition from Gaussian to power-law (~ 1/t 3 / 2 ) when 
the Heisenberg exchange coupling J becomes nonzero, 
and acquires a universal phase shift of 37r/4. The os- 
cillation frequency and phase shift as a function of time 
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can be used to determine the exchange and Overhauscr 
field fluctuations. We have also investigated the effects 
of an inhomogencous polarization on Ct (t) , and have 
suggested a general strategy for reducing decoherence in 
this system. Finally, we have discussed orbital dephasing 
and its effect on singlet-triplet decoherence. 
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APPENDIX A: ESTIMATING THE 
OVERHAUSER FIELD 

In this appendix we estimate the size of the Overhauscr 
field inhomogencity for a typical double quantum dot, 
and show that this quantity depends, in a sensitive way, 
on the form of the orbital wavefunctions. 

As in the main text, we take the average Overhauscr 
field and the Overhauser field inhomogeneity to be h = 
\ (hi + and 8h = | (hi — h2) respectively, where 

h; = AvJ^k |V'o( r fc) | 2 Ifc' an d ^o( r ) i s orbital eigenstate 
I in the double quantum dot. In the presence of tun- 
neling, the eigenstates of a symmetric double quantum 
dot will be well-describedS^ by the symmetric and an- 
tisymmetric linear combination of dot-localized states 
<Mr), 1 = 1,2: ^o' 2 (r) = 75 (<M r ) ± M 1 ))- In this case > 
we find 

(Sh) rms = Av(^Re[<Pl(v k )Mrk)]Ik) ■ (Al) 

\ k I rm s 

We take (i £ fc I fc > » y/I(I +1)/N to be the r.m.s. 
value for a system of N nuclear spins with uniform polar- 
ization p <C 1. Changing the sum to an integral according 
to v ^2 k — > / d 3 r then gives 

(6h) rms ~ 7] f^p±A = 1 a, (A2) 

where 7 = J d 3 r Re [0*(r)02(r)] is the overlap of the lo- 
calized orbital dot states and we have introduced the 
energy scale a = \J 1(1 + l)A/\HSf . The result in Eq. 
(|A2|) suggests that the Overhauser field inhomogeneity 
can be drastically reduced in a symmetric double quan- 
tum dot simply by separating the two dots, reducing the 
wavefunction overlap. If, however, the double dot is suffi- 
ciently asymmetric, the correct orbital eigenstates will be 



well-described by localized states ip l (r) = </>;(r), I = 1, 2, 
(with overlap 7 <C 1), in which case we find 

(Sh) imB ^ ] [^^-A = a. (A3) 

Thus, great care should be taken in determining (<5h) rms 
based on microscopic parameters. In particular, for a 
symmetric double quantum dot, the overlap 7 must also 
be known to determine (<5h) rms based on N. 

In contrast, for the total Overhauser operator h, in 
both of the above cases (ip ' (r) = 4^ (0i(r) ± 02 (r)) or 

^(r)=0,(r), 1 = 1,2), we find 

(h) rm8 «Y^^A = a. (A4) 

APPENDIX B: EFFECTIVE HAMILTONIANS 
FOR TWO-ELECTRON STATES IN A DOUBLE 
QUANTUM DOT 

In this appendix we derive effective Hamiltonians for 
a two-electron system interacting with nuclear spins in 
a double quantum dot via the contact hyperfine interac- 
tion. 

We begin from the two-electron Hamiltonian in second- 
quantized form, 

H = H S P + H c + H T + H z +H h{ , (Bl) 

where H$p describes the single-particle charging energy, 
He models the Coulomb interaction between electrons 
in the double dot, Ht describes tunneling between dot 
orbital states, Hz gives the electron Zeeman energy (we 
neglect the nuclear Zeeman energy, which is smaller by 
the ratio of nuclear to Bohr magneton: /^at/^b ~ 1CP 3 ) 
and Hhf describes the Fermi contact hyperfine interaction 
between electrons on the double dot and nuclei in the 
surrounding lattice. Explicitly, these terms are given by 

H S p = ^Vginiu; m < j=d\ a di a , (B2) 

la 

H c = U^n^mi + U'(n^ +n u )(n 2T +n 2i )(B3) 
1 

H T = tia {4ad2a + 4^i CT ) , (B4) 

(T 

Hz = |5> iT -n u ), (B5) 
1 

H M = J2 S ' ' h *; Sl = \j2 dia aa ,d la ,. (B6) 

l aa' 

Here, d\ creates an electron with spin a in orbital state 
I (I = 1, 2), V g i is the single-particle charging energy for 
orbital state I, U is the two-particle charging energy for 
two electrons in the same orbital state, and U' is the 
two-particle charging energy when there is one electron 
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in each orbital. When the orbital eigenstates are local- 
ized states in quantum dot I = 1,2, V g i is supplied by 
the back-gate voltage on dot I and U (W) is the on-site 
(nearest-neighbor) charging energy. t\2 is the hopping 
matrix element between the two orbital states, e z is the 
electron Zeeman splitting, h/ is the nuclear field (Over- 
hauser operator) for an electron in orbital I, and o a a' 
gives the matrix elements of the vector of Pauli matrices 
o = ("xj <Jy, a z ). In the subspace of two electrons occupy- 
ing two orbital states, the spectrum of Hsp + Hq consists 
of four degenerate "delocalized" states with one electron 
in each orbital, all with unperturbed energy £7(i,i) (a sin- 
glet 15(1, 1)) and three triplets: \Tj(l, 1)) ; j = ±,0), and 
two non-degenerate "localized" singlet states |S(2,0)) 
and |S(0, 2)), with two electrons in orbital I — 1 or / = 2, 
having energy £7(2, o) and E( 2 ), respectively. 

To derive an effective Hamiltonian H e g from a given 
Hamiltonian H, which has a set of nearly degenerate lev- 
els {\i)}, we use the standard procedure^, 

H cS = PHP + PHQ E _ 1 QHQ QHP, (B7) 

where P = J^. |i) (i\ is a projection operator onto the 
relevant subspace and Q = 1 — P is its complement. 
We choose the arbitrary zero of energy such that 
= V g i + V g 2 + U' = and introduce the detun- 
ing parameters 

Si = - £( 2 .o) = -2V g i -U = -S-U- C/(B8) 

Si = £(i,i) - £7 ( o, 2) = ~2V g2 -U = S. (B9) 

We then project onto the four-dimensional subspace 
formed by the delocalized singlet 15(1,1)) and three 
delocalized triplet states |Tj(l,l)), j = ±,0. That 
is, we choose Q = |5(0,2)) (5(0, 2)| + |5(2,0)) (5(2, 0)|, 
P = 1 - Q. When Si,S 2 > t i2 , we have E r* £7(i,i) = 
in the denominator of Eq. I|B7|) . This gives an effective 
spin Hamiltonian in the subspace of one electron in each 
orbital state: 

H eS = e 2 ^5f +^h,.S,- jQ -Si -S^BIO) 

J«-2t? 2 (i-- V (Bll) 

This Hamiltonian is more conveniently rewritten in terms 
of the sum and difference vectors of the electron spin and 
Ovcrhauscr operators S = Si + S2, SS = Si — S2 and 
h=i(hi+h 2 ), Sh= i(h x -h 2 ): 

H oS = e 2 S 2 + h-S + <5h-5S + ^S-S- J. (B12) 

Neglecting the constant term, in the ba- 
sis of singlet and three triplet states, 
{|S(1,1)) = |S) ,1^(1,1)) = \Tj) , j = ±,0}, the Hamil- 
tonian matrix for H e g takes the form 



/ -6h+/V2 Sh z 5h-/V2 \ 
-6h-/V2 J + e z +h z h~/V2 
Sh z h + /V2 J h-/V2 

\ 5H+/V2 h + /V2 J-e z -h z ) 

(B13) 

where Sh^ = Sh x ±i5h v and /i ± = h x ±ih y . We are inter- 
ested in this Hamiltonian in two limiting cases, where it 
becomes block-diagonal in a two-dimensional subspace. 

1. Effective Hamiltonian in the \S) — |To) subspace 

Projecting H onto the two-dimensional subspace 
spanned by |T ) -► \t z = +1) and |5) -> \t z = -1), we 
find 

H = N + iv • r, (B14) 

where r = (t x ,t v ,t z ) is a vector of Pauli matrices. The 
leading and first subleading corrections to Hq in powers 
of l/e z arc (£T = H^ + + ■■■, £f« = + v<*>): 



v r 


J 

~ 2' 
= J, 
= 2Sh z , 

- i<^" + 


] + [8h-,Sh+]), 


(B15) 

(B16) 
(B17) 
(B18) 


z(l) 

V 


= ^(["--" + 


} - [8h-,Sh+]), 


(B19) 


u 


= L(sh+h--\ 


-8h~h + ) . 


(B20) 



Here, N x = (N x ,N y x ,N^), v x = {v x ,v x ,v* x ), N x = 
N x ±iN x , and v x = v x ±iv x . For a typical unpolarized 
system, we estimate the size of all subleading corrections 
from their r.m.s. expectation values, taken with respect 
to an unpolarized nuclear state. This gives 

(ffo (1) ) CB21) 

\ / rms V £ z / 

where a is given by a = +T)A/VN (for a GaAs 

quantum dot containing N w 10 5 nuclear spins, 1/a ~ 
5 ns) . We therefore expect dynamics calculated under 

Hq ^ to be valid up to timescales on the order of e z /a 2 » 
1/a, when e z 3> a. 

2. Effective Hamiltonian in the \S) — \T+) subspace 

When the Zeeman energy of the \T + ) triplet 
state approximately compensates the exchange, 
max((h) rms ,(5h) rms ,|A|) « J (where A = e z + J), 
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wc find an effective Hamiltonian in the subspace 
\T + )^\t* = +1),\S)^\t* = -1); 



H+=N+ + iv+ • r, 



(B22) 



where the leading and sublcading corrections in powers 
of 1/ J are 



N 



(0) 



1 



(A + h z ). 



AO) 



,+(o) 



2 

A + h z , 
-V25h + , 



(B23) 

(B24) 
(B25) 



AT« = -i- ( (6h z ) 2 + hh-5h+ + \h-h+ )(B26) 



2 J 



,+(!) 



1 



1 



- I (<to*) z + -5h-Sh+ - 2 h ~ h ) > ( B27 ) 



= -V2 



J 



(B28) 



Once again, wc estimate the influence of the sublcading 
corrections from their r.m.s. value with respect to a nu- 
clear spin state of polarization p -C 1, giving 



H 



(i: 



= 



(B29) 



We therefore expect the dynamics under to be valid 
up to time scales on the order of t ~ J /a 2 3> 1/a for 
J > a. 



APPENDIX C: ASYMPTOTICS 

1. C"t (oo) for J > 2a , J < 2<r 

In the limit of J — > 0, we perform an asymptotic expan- 
sion of the integral in Eq. i|14|) by separating the prcfac- 
tor into a constant piece and an unnormalized Lorentzian 
of width J/2: 



C(x) ~ 1 ( 1 (J/2)2 
C{X) -2[ 1 { J,2f + : 



(CI) 



The Gaussian average over the constant term gives 1/2 
and when J/2 <C <Ja, the typical x contributing to the 
Lorentzian part of Eq. (|14fl is x < J/2 <C Co, so wc 
approximate exp(— ^x 2 /(Jq) « 1 in the integrand of this 
term. Integrating the Lorentzian then gives the result 
in Eq. 115fl for J -c 2cto. In the opposite limit of 
J 3> 2<toj the Lorentzian is slowly- varying with respect to 
the Gaussian, and the prefactor can be expanded within 
the integrand C{x) ~ 2x 2 /J 2 . Performing the remain- 
ing Gaussian integral gives the result in Eq. i|15|) f° r 
J > 2cr - 

2. C^oii) for f -» oo 

To evaluate the integral in Eq. (|13fl a t long times 
when J 7^ 0, we make the change of variables u = 

A 2 + (x/ao) 2 — A, A = J/2<7o, i = 2<7ot, which gives 



6%(i/2* ) =- -±= ^ du^±^ exp |-i (u 2 + 2uA) + i( u + A)t} , (C2) 

A = J/2ct , t = 2a t. (C3) 
I 



At long times, the major contributions to this integral 
come from a region near the lower limit, where u < l/t. 
For t > max(l/A, 1) (i.e. i > max(l/ J l/2cr )), wc 



approximate the integrand by its form for u <C max(A, 1), 
retaining the exponential term as a cutoff. This gives 



iXt 



C%(t/2a ) 



%A 



duy/^e~ [x - d)u = - 



iXl 



2VA (A - it) 



3/2 ' 



(C4) 



When t 3> A (i.e. t ^> J/Actq), we expand the denomina- 
tor of the above expression, which gives the result in Eq. 



G3- 
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